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ABSTRACT 

It has recently been proposed that ultrahigh energy (;> 10 19 eV) cosmic rays 
(UHECRs) are accelerated by the blast waves associated with GRBs. We calculate 
the observed synchrotron spectrum from protons and energetic leptons formed in the 
cascades initiated by photopion production, taking into account 77 attenuation at 
the source. Normalizing to the emission characteristics of GRB 970508, we predict 
~ 10 MeV - 100 GeV fluxes at a level which may have been observed with EGRET 
from bright GRBs, and could be detected with the proposed GLAST experiment 
or with ground-based air Cerenkov telescopes having thresholds <^ several hundred 
GeV. The temporal decay of the UHECR-induced high-energy 7-ray afterglows is 
significantly slower than that of the lower-energy burst and associated synchrotron 
self-Compton (SSC) radiation, which provides a direct way to test the hadronic origin 
of a high-energy GRB afterglow. Besides testing the UHECR origin hypothesis, the 
short wavelength emission and afterglows can be used to probe the level of the diffuse 
intergalactic infrared radiation field or constrain redshifts of GRB sources. 



Subject headings: cosmic rays — gamma rays: bursts — radiation mechanisms: 
nonthermal 



1. Introduction 



The cosmological origin of GRBs has very likely been confirmed as a result of observations 
of GRBs with the Beppo-SAX mission. The good imaging of X-ray emission from GRBs has led 
to optical counterpart identification of GRB 970228 (van Paradijs et al. 1997) and GRB 970508 
(Djorgovski et al. 1997). In the case of GRB 970508, observations of absorption lines in the 
spectrum of the optical counterpart place a lower limit on its redshift of z = 0.835 (Metzger et al. 
1997). 

The blast wave model for GRBs has met with considerable success in explaining the time 
dependence of the X-ray and optical afterglows (e.g., Meszaros, Rees, & Wijers 1997; Vietri 1997a; 
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Waxman 1997). The release of ss 10 52 ergs of energy in a small volume results in the formation 
of a relativistically expanding pair fireball which transforms most of the explosion energy into 
the kinetic energy of baryons in a relativistic blast wave. Reconversion of the kinetic energy into 
radiation occurs due to collisions between different shells ejected from the central source or when 
the blast wave decelerates as it sweeps up matter from the external medium. 

Shocks formed by these processes can accelerate protons to very high energies. If a significant 
fraction of the observed GRB power is transformed into UHECRs, then the measured energy 
density of these particles is in rough agreement with the power produced by GRBs, taking into 
account that UHECRs travel <J 100 Mpc before losing a significant fraction of their energy due to 
photopion production with the cosmic microwave background radiation (Zatsepin & Kuzmin 1966; 
Greisen 1966). On the basis of such arguments, Waxman (1995a, b) and Vietri (1995) proposed 
that UHECRs are accelerated by GRBs. Associated ~ 10 14 eV neutrino production (Waxman &; 
Bahcall 1997) and GeV photon production (Vietri 1997b) due to UHECR acceleration during the 
prompt 7-ray emission phase of the GRB have been recently predicted. Detection of neutrinos in 
temporal and spatial association with GRBs would confirm this cosmic ray origin theory, but will 
require larger neutrino detectors than presently exist. 

In this Letter, we calculate both the prompt and delayed high energy 7-ray emission due to 
UHECRs accelerated in GRBs. We calculate the emergent synchrotron radiation from protons, 
from positrons produced in the decay of 7r + , and from pairs produced by 77 interactions. The 
observed emission characteristics of GRB 970508 are used to set parameter values for the blast 
wave model of GRBs (see Meszaros & Rees 1993; Meszaros, Rees Sz Papathanassiou 1994; Dermer 



After reaching a maximum speed defined by the baryon content and total energy of the 
explosion, the Lorentz factor V of a relativistic blast wave evolves by sweeping up matter from 
its surrounding medium (Blandford & McKee 1976; Meszaros & Rees 1993). For a blast wave of 
total energy E = W 52 E§2 erg, which expands with initial bulk Lorentz factor Tq = 300r,3oo into a 
medium of uniform density uq cm" 3 , the coasting phase ends and the deceleration phase begins 
when the blast wave has a characteristic radius given by 



(Rees & Meszaros 1992). For x > xo, the bulk Lorentz factor follows the power-law behavior 
T(x) = Tq {x/xq)~ 9 where, for a uniform density medium, g = 3/2 in the non-radiative regime 
when transformation of the internal energy of swept-up particles into radiation is inefficient, and 
g = 3 in the radiative regime where the internal energy is efficiently radiated. 



k Chiang 1998). 



2. 



Evolution of GRB Blast Wave 




(1) 
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According to the scenario of UHECR acceleration by GRB blast waves (Vietri 1995; Waxman 
1995a,b), we assume that a fraction £ of the energy in nonthermal protons in the blast wave region 
is transformed into a power-law distribution iV cr (7 cr ;x) = K cr (x)'y~ s of ultrarelativistic cosmic 
ray protons with Lorentz factors T(x) < 7 cr < 7 cr ,max in the fluid frame comoving with a small 
element of the blast wave region. Here 



K (x) - — ^ pK ' v -1— (2) 



1 jE p (x){2- s) 
m p c 2 7^7^ - T(x) 

where E p (x) is the energy in non-thermal protons swept up from the external medium, and the 
maximum cosmic-ray Lorentz factor is 



^max-lO^r-^V/V/ 2 (3) 

(Vietri 1995, 1997b). The term ( in this equation is the equipartition factor for the magnetic field, 
given through the expression 

H(G) = (87rrm p c 2 n 

)i/2 r ^i/2 ^ a6 r[( r /i ) n C] 1/2 , (4) 

and r is the compression ratio which can exceed 4 for relativistic shocks or in shocks where 
nonlinear feedback of the nonthermal particle pressure affects the shock structure (e.g., Ellison, 
Jones, & Reynolds 1990). We let £ = 1, r = 10, and £ = 1 in our calculations. The choice of 
( ~ 1 is needed for Fermi acceleration to be efficient enough to produce cosmic-ray protons of 
E p ~ 10 20 eV (Vietri 1995). 

Finally, we have to relate the measured flux at a given time to the comoving differential 
photon density n p h(e), where e = Ej/(m e c 2 ) is the dimensionless photon energy. The observed 
radiation is a convolution of the contributions from regions of the blast wave emitting at 
different times and moving with different Lorentz factors and directions with respect to the 
observer (Rees 1967). The completely self-consistent transformation to the observer's frame 
can therefore in general only be done numerically, taking into account the evolution of the 
intrinsic radiation spectrum. However, due to the very strong Doppler enhancement along the 
direction of motion in a relativistic blast wave, only a small fraction of the blast wave covering 
a solid angle of order 1/r 2 contributes significantly to the observed radiation, within which 
the difference in light travel time to the observer is negligible. Within this approximation, an 
observed power-law spectrum given by S'obs( € obs) = So e~ b " between e b s ,i < e bs < e obs,2 is related 
to the spectral photon density in the comoving frame, given by n p h(e) = 7ip h e~( 1+a ) between 
e 1 = (1 + z) e obs ,i/(2r) < e < (1 + z) e obSi2 /(2r) = e 2 , through the expression 



ph ~ cx 2 m e c 2 (2T) 1 + a ' { ) 



Here, di is the luminosity distance. 
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3. Photopion Production, Proton Synchrotron Radiation and Pair Cascades 

We calculate the 7-ray and positron production spectra from decaying pions produced by 
UHECRs interacting with photons (see, e.g., Stecker 1979). The UHECR energy- loss rate through 
photopair production is much smaller than that of photopion production for the proton energies 
and photon spectra considered here, and will be neglected. The differential cross section for 
photopion production is approximated by da/de' ~ <To<5(e' — E A /m e c 2 ), where o"o = 2 • 10 -28 cm 2 , 
e' = 7 cr e(l — n) is the photon energy in the proton rest frame, and E A = 330 MeV. Noting that 
pions are produced primarily near threshold in the protons' rest frame (i. e., j n ps 7 cr ), we obtain 
the pion production rate 
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The generalized Heaviside function B is defined as @(x;a, b) = 1 if a < x < b, and otherwise. 
The photon spectrum from decaying neutral pions is then 
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For simplicity, we assume that charged pions are produced at the same rate (Eq. [6]) as ir°s 
and, after decaying into one positron and three neutrinos, on average 1/4 of the energy of the 
pion is carried by the positron (instead of 1/2 carried by each of the two photons in the case of 
7T° decay). This yields iV e +(7 + ) iV^ ^ 27 (27+) for the first generation of relativistic positrons 
injected into the shock region. We take into account energy losses due to synchrotron radiation, 
but neglect escape, pair annihilation, and Compton losses. Escape is unlikely given the small 
Larmor radius, and pair annihilation is only important for low-energy positrons which radiate 
inefficiently. Compton losses are found to be much less important than synchrotron losses for the 
parameters used here. Integration of the continuity equation yields the "steady-state" distribution 
of the first generation of positrons. 

The proton synchrotron spectrum is calculated using a (5-function approximation for the 
emission of a single proton, giving 
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where e H , P = (3/2)(m e /m p ) H/H c , and F c = 4.414 • 10 13 G. 

We compute the opacity r 77 using the full pair production cross section and assuming that the 
shell width (5x s h = x/T (Meszaros, Laguna, & Rees 1993). The shock region is optically thick to 
77 pair production on the soft burst and afterglow radiation for photons with Ey ^ 100 GeV. This 
differs from the results of Waxman & Bahcall (1997) and Vietri (1997b), who find the region to be 
optically thick for photons above ~ 100 MeV and ~ 1 TeV, respectively. This difference is largely 
due to the smaller size scale used by Waxman & Bahcall and the different nominal values of the 
Lorentz factors (I^oo = 1/3, 1, and 3 in Waxman h Bahcall, this study, and Vietri, respectively). 
The 77 absorbed portion of the high-energy photons initiates a pair cascade. Since most of the 
secondary pairs of the electromagnetic cascade are injected at high energies, the steady-state 
pair distribution of secondary particles can be approximated by a power-law with index 2 over 
a wide range of electron/positron energies. This yields a synchrotron spectrum of spectral index 
a syn = 1/2. We approximate the radiative output of the cascade by a 77-absorbed power-law 
spectrum of energy index 1 /2 containing the energy of the absorbed high-energy photons from 7r° 
decay and synchrotron radiation of the first generation of positrons from ir + decay. 

SSC radiation from the primary electrons may also play an important role. The numerical 
simulations described in Chiang & Dermer (1998) were used to determine the expected level and 
energy range of the SSC component. Note that the SSC component is significantly weaker than 
estimated using simple Thomson-limit arguments because Klein-Nishina effects are very important 
in this case. 



For the purpose of illustration, we compute the expected high-energy prompt radiation and 
afterglow from GRB 970508. We approximate its observed spectrum during the burst phase by a 
broken power-law with energy index a± = below inbreak ~ 45 keV and a<i = 1 above the break, 
fitting to the Beppo-SAX fluxes measured in the 2-26 keV and 40-700 keV bands during the first 
25 s of the GRB (Piro et al. 1998). The spectrum is assumed to extend from E^bs,! = 1 eV, 
corresponding to the self-absorption frequency early in the GRB, to -Eobs,2 = 10 MeV (the 
calculations do not depend sensitively on the upper energy). We assume E = 10 52 erg, z = 0.835, 
To = 300. Eq. (3) implies that the cosmic-ray spectrum extends up to 7 C r,max = 7 • 10 8 , and we 
assume that it has a spectral index of s = 2. Fig. 1 shows the different components of the resulting 
broad-band spectrum at the point xo, and the pair production opacity r 77 at this point. 

Examing synchrotron losses in the evolving magnetic field (eq. [4]), one finds that only 
particles with 7 <; 5(m/m e ) 3 /[(r/10)Cno^ 3 r3QQi?52 3 ] cool efficiently through synchrotron proceses. 



4. 



Numerical results 



-6- 



Thus synchrotron cooling is very inefficient for protons, though electrons cool efficiently if ( ss 1. 
The radiative output of protons is dominated by synchrotron radiation rather than photo-pion 
production, although the efficiency of the latter process can be increased in a colliding shell 
scenario (e.g., Waxman & Bahcall 1997). 

Approximating the observed evolution of the optical to soft 7-ray spectrum of the GRB and 
the afterglow from GRB 970508, we assume that the flux remains constant for to = 30 s in the 
observer's frame and then decays as S(t) = So {t/to)~ im ■ The break energy -Screak and the upper 
cut-off E2 are assumed to decrease with the same temporal behavior. Because cooling is inefficient, 
we expect that the evolution of the blast wave is more closely described by the non-radiative 
limit. Fig. 2 illustrates the time dependence of the high-energy 7-ray afterglow, compared to the 
evolution of the low-energy burst and the SSC radiation in this case. The high-energy afterglow 
decays much more slowly than the optical-to-X-ray afterglow and the SSC component. In 
the radiative case, the proton synchrotron afterglow radiation decays much faster than in the 
non-radiative case, but still considerably less rapidly than the optical-to-X-ray afterglow. 

The temporal behavior of the uncooled proton synchrotron flux in the evolving field (4) goes 
as i~ b * where Xuncooicd = [4<7 + (7r//4) — 3]/ (1 + 2g) for s = 2 and an external medium n ext oc x~ v 
(see Dermer & Chiang 1998, eqs. [50], [51]). In the cooled regime, characteristic of the electon 
synchrotron radiation, we find Xcooled = (4g + 77 — 2)/(l + 2g). In the non-radiative g = 3/2 regime 
with 77 = 0, WG find that Xuncooled — 

3/4 and 

Xcooied — 1- The effects of the cascade cause the 
high-energy radiation to decay even more slowly than estimated by these analytic expressions. 

The level of the proton synchrotron flux is not strongly dependent on the endpoints of the 
low-energy spectrum because this component is produced mainly at energies where intrinsic 77 
attenuation is unimportant. The photopion component is, however, reduced when the high-energy 
cutoff of the burst emission is decreased. Due to the smaller 77 optical depth, a smaller portion 
of the high-energy radiation resulting from tt° decay and synchrotron emission of first-generation 
positrons is absorbed and transferred to the cascade. Different values of To significantly change 
the energy range and amplitude of the emergent proton synchrotron and photopion emission, but 
only during the early portion of the GRB. At comparable times later in the GRB, the flux and 
energy ranges are similar. 

5. Discussion 

The high energy 7 radiation from proton synchrotron and the photopion-induced pair cascade 
decays more slowly than the low energy afterglow and the SSC radiation, because cooling is 
inefficient for protons whereas electrons cool efficiently. Therefore, synchrotron radiation from 
protons can be distinguished from direct electron synchrotron or SSC emission by comparing the 
temporal decay rates in the afterglow phase, even though the SSC radiation is produced at roughly 
the same flux levels as the proton synchrotron radiation early in the GRB. A monotonically 
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increasing ratio of the high-energy 7-ray flux to the optical or soft X-ray flux in the GRB afterglow 
would provide strong evidence for proton synchrotron radiation and UHECR acceleration in 
GRBs. 

The flux sensitivities of EGRET and GLAST (pointed mode) for an a = 1 spectrum are 
« 8 • 1CT 8 ph (> 100 MeV) cm" 2 s" 1 and « 2 • 10~ 9 ph (> 100 MeV) cm' 2 s _1 , respectively 
(Kurfess et al. 1997; GLAST 1997). In the inset to Fig. 2, we compare the fluence sensitivities of 
EGRET, GLAST and Whipple to the fluence predicted by our calculations of GRB 970508. The 
calculated > 100 MeV burst fluence is an order of magnitude below the EGRET sensitivity level, 
but could be detected in brighter GRBs. GLAST, when operated in the pointed mode, would be 
able to detect the > 100 MeV 7-ray afterglow at the predicted level. The emission from UHECRs 
accelerated in GRBs may already have been observed in six bright BATSE bursts that were 
detected with EGRET (Dingus 1995; Catelli, Dingus, & Schneid 1996; Hurley et al. 1994) and 
display weak evidence for a hard tail at > 100 MeV and GeV energies. From Fig. 2, we note that 
the best opportunity for detection of ~ 100 MeV emission with EGRET or GLAST is within the 
first several minutes of the GRB whereas a larger fraction of the high-energy gamma-ray emission 
is produced on time scales of hours, which is in accord with the detection of an 18 GeV photon 
from GRB 940217 90 minutes after the main part of the GRB (Hurley et al. 1994). 

Fig. 2 indicates weak delayed > 100 GeV-TeV emission, which is below the thresholds of 
existing Gerenkov telescopes. This emission would not, in any case, be detectable from a GRB 
located at a redshift z ^ 0.5 (see Stecker &; de Jager 1997; Salamon & Stecker 1998) due to 
absorption by the intergalactic infrared radiation field, but could be detected from nearby bright 
GRBs. The delayed afterglows in the 10-100 GeV range from such GRBs could provide a way 
to measure the level of the infrared background radiation and to yield upper limits on redshifts, 
but will require high sensitivity ground-based Gerenkov telescopes with low energy thresholds. If 
redshifts are independently measured, cutoffs in the observed high energy radiation spectrum can 
be used to infer the level of the infrared background. 

The photopion process could be more important during the early phase of the GRB if the 
size scale of the system is determined by collisions between shells (Waxman & Bahcall 1997). In 
this case, we would expect a considerably larger > 10 MeV flux due to the photopion/77 cascade 
during the primary 7-ray emission phase. Flux levels of high-energy 7 radiation could therefore 
constrain the expected level of neutrino production from UHECRs and discriminate between the 
colliding shell and external shock scenario of GRBs (Vietri 1997b). 
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Fig. 1. — Broadband spectrum of a GRB blast wave at the point xq of transition to the asymptotic 
regime, for parameters derived from the observations of GRB 970508. Dashed: optical to soft 
7-ray spectrum; dotted: proton synchrotron radiation; long dashed: synchrotron radiation from 
positrons produced in ir + decay; dot-dashed: ir° decay; solid: total high-energy spectrum, including 
77 absorption and synchrotron emission from the pair cascade; thick dashed: 77 pair production 
opacity multiplied by 10~ 10 . 
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Fig. 2. — Time history of the broadband burst and afterglow emission in the non-radiative regime 
at different observer times. The expected SSC radiation associated with the low-energy burst 
emission is shown for comparison. The inset shows time- and energy-integrated fluxes for high 
energy emission from GRB 970508 for energies > 100 MeV (dotted curve) and > 350 GeV (short- 
dashed curve). The dot-dashed, solid and long-dashed line are the fluence sensitivities for EGRET, 
GLAST in its pointed mode, and Whipple, respectively. 



